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We present a Monte Carlo simulation of the perturbative Quantum Chromodynamics (pQCD) 
shower developing after a hard process embedded in a heavy-ion collision. The main assumption is 
that the cascade of branching partons traverses a medium which (consistent with standard radiative 
energy loss pictures) is characterized by a local transport coefficient q which measures the virtuality 
per unit length transferred to a parton which propagates in this medium. This increase in parton 
virtuality alters the development of the shower and in essence leads to extra induced radiation and 
hence a softening of the momentum distribution in the shower. After hadronization, this leads 
to the concept of a medium-modified fragmentation function. On the level of observables, this is 
manifest as the suppression of high transverse momentum (Pr) hadron spectra. We simulate the soft 
medium created in heavy-ion collisions by a 3-d hydrodynamical evolution and average the medium- 
modified fragmentation function over this evolution in order to compare with data on single inclusive 
hadron suppression and extract the q which characterizes the medium. Finally, we discuss possible 
uncertainties of the model formulation and argue that the data in a soft momentum show evidence of 
qualitatively different physics which presumably cannot be described by a medium-modified parton 
shower. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 

Jet quenching, i.e. the energy loss of hard partons cre- 
ated in the first moments of a heavy ion collision due to 
interactions with the surrounding soft medium has long 
been regarded a promising tool to study properties of 
the soft medium [J, 0, H, |j, [E El • The basic idea is to 
study the changes induced by the medium to a hard pro- 
cess which is well-known from p-p collisions. A number 
of observables is available for this purpose, among them 
suppression in single inclusive hard hadron spectra Raa 
Q , the suppression of back- to-back correlations [1, or 
single hadron suppression as a function of the emission 
angle with the reaction plane • 

Calculations have now reached a high degree of sophisti- 
cation. Different energy loss formalisms are used together 
with a 3-d hydrodynamical description of the medium 
[Tl| in central and noncentral collisions to determine the 
pathlength dependence of energy loss [HI, [IH, 0, fill ]. 
Some of these models have also been employed success- 
fully to describe the suppression of hard back-to-back 
hadron correlations [ill [l?], EH • 

The existing formulations of energy loss can roughly be 
divided into two groups: Some compute the energy loss 
from a leading parton 0, [H, Eil whereas others com- 
pute an in-medium fragmentation function by following 
the evolution of a parton shower [l3l [20j . Recently, the 
Monte Carlo (MC) code JEWEL [2l| has also been de- 
veloped which simulates the evolution of a parton shower 
in the medium in a non-analytic way. This model builds 
on the success of MC shower generators like PYTHIA 
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[H [H or HERWIG for showers in vacuum. 

In the present work, we follow an approach which is very 
similar to the one taken with JEWEL, i.e. we modify a 
MC code for vacuum shower to account for medium ef- 
fects. However, while JEWEL so far chiefly implements 
elastic scattering with medium constituents and includes 
radiative energy loss only in a schematic way, we rather 
wish to focus on radiative energy loss in the following. 
This is based on the observation that elastic energy loss 
has the wrong pathlength dependence to account prop- 
erly for the suppression of back-to-back correlations [25| 
and hence cannot be a large contribution to the total en- 
ergy loss of light quarks. In particular, we assume that 
partons traversing the medium pick up additional virtu- 
ality which induces additional branchings in the shower, 
thus softening the parton spectrum, but that there is no 
transfer of longitudinal momentum from the hard parton 
to the scatterers in the medium. 

This work is organized as follows: First, we outline the 
computation of hadron spectra in the formalism, starting 
from the hard process. The key ingredient of our model, 
the medium-modified fragmentation function (MMFF) is 
described in detail in section IIIII where we outline the 
MC simulation of showers in vacuum and present how 
the algorithm is modified to simulate showers in medium. 
We present various observables which show the expected 
modification of the jet by the medium. In section fVl we 
present a comparison of the suppression calculated us- 
ing the MMFF in a 3-d hydrodynamical model for the 
medium evolution ,11} with the measured nuclear sup- 
pression in central AuAu collisions at 200 AGeV and use 
this result to extract an estimate for the medium trans- 
port coefficient q. We follow with a discussion of the 
model uncertainties. Finally the limits of the approach 
in the light of the patterns seen in semi-hard and soft 
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correlations with a hard trigger hadron are discussed. 

II. THE HARD PROCESS 

We aim at a description of the production of high Pt 
hadrons both in p-p and in Au-Au collisions. The un- 
derlying hard process can be computed in leading order 
(LO) pQCD. We assume in the following that the ac- 
tual hard process is not influenced by the fact that a 
soft medium is created in Au-Au collisions, but that the 
subsequent parton shower (which extends to timescales 
at which a medium is relevant) is modified by the pres- 
ence of such a medium, whereas hadronization itself takes 
place sufficiently away from the medium such that it can 
be assumed to take place as in vacuum. In this section, 
we describe the computation of the hard process itself. 
The production of two hard back to back partons k, I with 
momentum px in a p-p or A-A collision in LO pQCD is 
described by 



d(J AB^kl+X ^frij^kl 
, 2 , , = y2x 1 f i/A (x 1 ,Q 2 )x 2 f : j/B(x2,Q 2 ) £ 

dp T dyxdy 2 ~" ' dt 

(1) 

where A and B stand for the colliding objects (protons 
or nuclei) and yu 2 ) is the rapidity of parton k{l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction X\ and a factorization scale Q ~ px is 
fi/A( x i>Q )■ The distribution functions are different for 
the free protons [IE H3| and protons in nuclei [H, [2^ . 
The fractional momenta of the colliding partons i, j are 
given by xi >2 = ^ (exp[±j/i] + exp[±y 2 ])- 

Expressions for the pQCD subprocesses — (s, i, u) as 
a function of the parton Mandelstam variables s, t and u 
can be found e.g. in 30]. Inclusive production of a parton 
flavour / at rapidity y/ is found by integrating over either 
2/i or 2/2 and summing over appropriate combinations of 
partons, 



da AB^f+X 

dp T dy f 



1 1 



(ij),{kl) 



kl J- "I" Oij 



xif l /A(xi,Q 2 )x2f j /B(x 2 ,Q 2 ) 



+xifj/A(xi,Q 2 )x 2 f i/B (x 2 ,Q 2 ) 



(s,t,u)5fk H -7 — (s,u,t)5fi 



dt 1 dt 

— — — (s,u,t)8 fk + — — — (s,t,u)6fi 
dt dt 



(2) 



where the summation (ij),(kl) runs over pairs 
gg, gq, gq^qq,qq,7jq and q stands for any of the quark 
flavours u, d, s. 

For the production of a hadron h with mass M , trans- 
verse momentum Pt at rapidity y and transverse mass 
ttit — y M 2 + P T from the parton /, let us introduce 
the fraction z of the parton energy carried by the hadron 
after fragmentation with z = Eh/ Ef. Assuming collinear 
fragmentation, the hadronic variables can be written in 
terms of the partonic ones as 



rtiT cosh y — zpx cosh y f and ttit sinh y = Pt sinh y f . 

(3) 

Thus, the hadronic momentum spectrum arises from 
the partonic one by folding with the distribution 
D (z, [i 2 ) which describes the average number of frag- 
ment hadrons carrying a fraction z of the parton momen- 
tum at a scale fif ~ Pt as 



J 



da AB^h+X r d(J AB^f+X A 

= 2^ / dp 2 T dy f ^ 2 j dzD f ^ h {z,p 2 f )8 (m 2 T - M T {p T ,y f ,z)) 8(y -Y(jp T ,y { ,zj) (4) 



dP 2 dy 



with 



and 



/ P T 

Y(pt, 2/f , z) = arsinh sinhi/f 

V tut 



(6) 



M T (p T , yj, z) = (zp T ) 2 + M 2 tanh 2 y f (5) The lower cutoff z min = ^f- cosh?/ arises from the fact 
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that there is a kinematical limit on the parton momen- 
tum; it cannot exceed y/s/(2 coshy/) and thus for given 
hadron momentum there is a minimal z corresponding to 
fragmentation of a parton with maximal momentum. 
In the following, we describe this key ingredient in our 
formulation of energy loss, i.e. the computation the frag- 
mentation function both in vacuum and in medium using 
a pQCD shower evolution code. 



III. THE MEDIUM-MODIFIED 
FRAGMENTATION FUNCTION 

In this section, we describe how the medium-modified 
fragmentation function is obtained from a computation 
of an in-medium shower followed by hadronization. Key 
ingredient for this computation is a pQCD MC shower 
algorithm. In this work, we employ a modification of the 
PYTHIA shower algorithm PYSHOW In the ab- 

sence of any medium effects, our algorithm therefore cor- 
responds directly to the PYTHIA shower. Furthermore, 
the subsequent hadronization of the shower is assumed 
to take place outside of the medium, even if the shower 
itself was medium-modified. It is computed using the 
Lund string fragmentation scheme [3l| which is also part 
of PYTHIA. 



A. Shower evolution in vacuum 

We model the evolution from some initial, highly vir- 
tual parton to a final state parton shower as a series of 
branching processes a — > b + c where a is called the par- 
ent parton and b and c are referred to as daughters. In a 
longer chain of branchings, all partons which have been 
part of the evolution branch towards some parton i are 
called the ancestors of i in the following. 
In QCD, the allowed branching processes are q — > qg, 
g — > gg and g — » qq. The kinematics of a branching is 



described in terms of the virtuality scale Q 2 and of the 
energy fraction z, where the energy of daughter b is given 
by Eb = zE a and of the daughter c by E c — (1 — z)E a . It 
is convenient to introduce t — InQ 2 / Aqcd where Aqcd 
is the scale parameter of QCD. t takes a role similar to a 
time in the evolution equations, as it describes the evo- 
lution from some high initial virtuality Qo (to) to a lower 
virtuality Q m (i m ) at which the next branching occurs. 
In terms of the two variables, the differential probability 
dP a for a parton a to branch is [H, [33| 



dP a = ^2 ^P a ^ bc (z)dtdz 



b.c 



(7) 



where a s is the strong coupling and the splitting kernels 
Pa-.bc (z) read 



,(z)=4/3 



1 



P 



g^gg 



(z) = 3 



1 - z 
(l-z(l-z)) 2 
z{l-z) 



P g ^ q - q (z) = N F /2(z 2 + (1 - zf) 



(8) 

(9) 
(10) 



where we do not consider electromagnetic branchings in- 
volving a photon or a lepton pair and Np counts the 
number of active quark flavours for given virtuality. 
At a given value of the scale t, the differential probability 
for a branching to occur is given by the integral over all 
allowed values of z in the branching kernel as 



Ia^bc{t)= f Z+t ] dz^P a ^ bc ( Z ). (11) 
Jz-(t) 2lr 



The kinematically allowed range of z is given by 



J 



z± 




pj y/(MZ - Ml - M, 



AMjMl 



M 2 



(12) 



r 



where M 2 = Q 2 + mf with m< the bare quarm mass or 
zero in the case of a gluon. Note that the gluon emis- 
sion from a quark Eq. ((SJ) is singular at z = 1 and the 
gluon splitting probability ^ at z = 1 and z = 0, hence 
they lead preferably to soft gluon radiation close to the 
kinematic cutoff z± . Note also that a high virtuality Q a 
shrinks the available range of z and hence enforces on 
average harder radiation. 

Given the initial parent virtuality Q 2 or equivalently t a , 
the virtuality at which the next branching occurs can 



be determined with the help of the Sudakov form fac- 
tor S a (t), i.e. the probability that no branching occurs 
between to and t m , where 



S a (t) = exp 



Co 



dt'J2la^bc(t') 



b.c 



(13) 



Thus, the probability density that a branching of a occurs 
at t m is given by 
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dPg 

dt 



b.c 



S a (t). 



(14) 



These equations are solved for each branching by the 
PYSHOW algorithm [23| iteratively to generate a shower. 
For each branching first Eq. (fT4|) is solved to determine 
the scale of the next branching, then Eqs. ([8|)- (fT0|) are 
evaluated to determine the type of branching and the 
value of z, if the value of z is outside the kinematic 
bound given by Eq. (fP2|) then the event is rejected. Given 
to,t m and z, energy-momentum conservation determines 
the rest of the kinematics except for a radial angle by 
which the plane spanned by the vectors of the daughter 
parents can be rotated. 

In order to account in a schematic way for higher order 
interference terms, angular ordering is enforced onto the 
shower, i.e. opening angles spanned between daughter 
pairs b, c from a parent a are enforced to decrease ac- 
cording to the condition 



z b (l - z)b) > 1_ 



M 2 



M 2 



(15) 



After a branching process has been computed, the same 
algorithm is applied to the two daughter partons treating 
them as new mothers. The branching is continued down 
to a scale Q m in which is set to 1 GeV in the MC simula- 
tion, after which the partons are set on-shell, adjusting 
transverse momentum to ensure energy-momentum con- 
servation. 



B. Shower evolution in medium 

We now embed this shower evolution in a soft medium. 
It is important to realize that this implies that energy 
and momentum in this case is in general not conserved if 
one includs only one component (shower or medium) in 
the simulation, as there is flow of energy and momentum 
through interactions from the shower to the medium and 
vice versa. 

In this particular model, we aim at a description of radia- 
tive energy loss, i.e. we assume a medium which does not 
absorb momentum by recoil of its constituents, but rather 
increases the virtuality of partons propagating through 
it and thus inducing additional radiation. While this 
is clearly an idealized assumption of the medium, it al- 
lows to study cleanly the effect of induced radiation on 
a shower, and moreover it is expected to capture the 
physics of the dominant component of energy loss [25| . 
Such a medium can be characterize by a transport coef- 
ficient q which represents the increase in virtuality AQ 2 
per unit pathlength of a parton traversing the medium. 



Note that this represents an average transfer, i.e. a pic- 
ture which would be realized in a medium which is char- 
acterized by multiple soft scatterings with the hard par- 
ton. Alternatively, a model in which hard scatterings 
which occur probabilistically transfer the virtuality could 
be implemented. 

While the vacuum shower develops only in momentum 
space, the interaction with the medium requires mod- 
elling in position space as well. In order to make the 
link, we assume that the formation time of a shower par- 
ton with virtuality Q is developed on the timescale 1/Q, 
i.e. the lifetime of a virtual parton with virtuality Q 
coming from a parent parton with virtuality Q a is in the 
rest frame of the original hard collision (the rest frame 
of the medium may be different by a flow boost as the 
medium may not be static) given by 



n 



Eb_Eb_ 

Qt Ql 



(16) 



The time at which a parton a is produced in a branch- 
ing can be determine by summing the lifetimes of all an- 
cestors. Thus, during its lifetime, the parton virtuality 
is increased by the amount 



T +T„ 



dtq(Q 



(17) 



where £ is integrated along the spacetime path of the 
parton through the medium and q(£) = q{r, r, <p, rj s ) de- 
scribes the spacetime variation of the transport coeffi- 
cient where q is specified at each set of coordinates proper 
time r, radius r, angle <p and spacetime rapidity n s . In 
the following, we will use a hydrodynamical evolution 
model of the medium [llT | for this dependence. If the par- 
ton is a gluon, the virtuality transfer from the medium 
is increased by the ratio of their Casimir color factors, 
3/f = 2.25. 

If AQ 2 -C Q 2 , holds, i.e. the virtuality picked up from 
the medium is a correction to the initial parton virtuality, 
we may add AQ 2 to the virtuality of parton a before 
using Eq. (fl4]) to determine the kinematics of the next 
branching. If the condition is not fulfilled, the lifetime 
is determined by Q 2 + AQ 2 and may be significantly 
shortened by virtuality picked up from the medium. In 
this case we iterate Eqs. (TTB1) . (JTT|) to determine a self- 
consistent pair of (r a ,AQ 2 ). This ensures that on the 
level of averages, the lifetime is treated consistently with 
the virtuality picked up from the medium. 
The need for the self-consistent iteration must be investi- 
gated a posteriori when an actual q(() is specified. At this 
point, we note that it may frequently occur in a constant 
medium that virtuality transferred from the medium de- 
termines the lifetime of a parton, as towards the end of 
the shower development Q 2 may become small, leading to 
large lifetimes and large AQ 2 . But the situation is much 
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less severe for an expanding medium in which q ~ l/r" 
where a > 1 and r is the lifetime of the medium. 
We note that in this model, the shower gains energy from 
the medium. This is not unexpected, as we have only 
included energy transfer to the shower, but no interac- 
tion which would feed energy from the shower into the 
medium. In principle, one may argue that sufficiently 
soft partons of the shower would thermalize and thus be- 
come part of the medium, but at the moment we refrain 
from modelling energy transfer to the medium. Instead, 
we note that the increase in parton virtuality leads to 
induced radiation, transferring energy from hard partons 
to soft radiation and thus to a softening of the momen- 
tum distribution of partons inside the shower, along with 
an increase in momentum \qj<\ transverse to the shower 
initiating parton direction. 

C. Properties of the hydrodynamical medium 

After using the Lund string fragmentation scheme j3lj to 
hadronize the shower, we can compute the longitudinal 
momentum distribution of a shower initiated by a par- 
ton at fixed energy and hence the fragmentation function 
into any hadron species at this given partonic momentum 
scale. 

It is clear that the form of the medium-modified frag- 
mentation function (MMFF) is dependent on q(Q and 
that no general form can be given without specifying this 
quantity. However, we are only interested in very specific 
forms of q(C), he. those which occur in the soft medium 
created in an heavy-ion collision. 
If we assume the relation 




T [fm/c] 



FIG. 1: Transport coefficient q encountered by a quark trav- 
elling into +x direction through the hydrodynamical medium 
in a 200 AGeV collision at impact paramter b = 2.4 fm for 
different initial hard vertex position in the transverse (x, y) 
plane (in all cases y = 0). Values have been extracted for 
K = 1.5 (see Eq. (|18[1 . which represents the best fit to the 
data. 



We thus investigate in the following three different sce- 
narios (approximately representing a parton travelling 
into +x direction originating from x = 4 fm (A), x = 
(B) and x — —4 fm (C), y = in all cases in the trans- 
verse (x, y) plane at midrapidity. These trajectories are 
characterized by the parameters (b — 1.5, c = 3.3, Tp = 
5.8 fm/c) (A), (b = 1.5, c= 2.2, t e = 10 fm/c) (B) and 
(b = 1.5, c = 2.2, te = 15 fm/c) (C) and are quite typical 
for partons close to the surface (A), emerging from the 
central region (B) or traversing the whole medium (C). 



q(Q = K- 2- [e(C)] 3/4 (cosh KC) - sinhp(C) cos^) (18) 



D. The medium-modified fragmentation function 



between q and the medium energy density e, the local 
flow rapid ity p w ith angle ip between flow and parton 
trajectory [3J,|35|, we find that the vast majority of paths 
leads to a g(C) which can be described by (see Fig. Q] for 
examples) 



<z(C) 



(b + T/(lfm/c)Y 



(19) 



The exception are paths close to the surface for which 
the parton travels inward. However, such events are not 
only suppressed by the overlap function which is small 
at the medium surface but also very suppressed due to 
the strong medium effect for such long paths. Thus, the 
medium which will be probed by a parton emerging from 
a particular vertex can be characterized by the four pa- 
rameters (a, b, c, te) where te is the time at which the 
parton emerges from the medium. In practice, the range 
of parameters can be sufficiently narrowed down to 1 
< b < 2, te < te (with Tp the freeze-out time of the 
medium) and 2 < c < 4. 



In Fig. [2] we present the longitudinal momentum q z dis- 
tribution of hadrons in a jet originating from a 20 GeV 
quark, both in vacuum and medium-modified according 
to the three scenarios outlined above. Instead of specify- 
ing the parameter a, we characterize the scenarios instead 
by the total amount of virtuality 



aqL = 



(20) 



which would be transferred to a single parton crossing 
the medium (note that since branchings occur, multiple 
partons travel through the medium in a shower, hence the 
virtuality transfer to the shower as a whole is larger than 
AQj ot , but can only be computed on an event-by-event 
basis) . 

In general, the medium modification induces the ex- 
pected changes to the shower, i.e. the high q z tail is 
reduced by the medium, whereas induced radiation cre- 
ates additional multiplicity at low q z . This transfer of 
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FIG. 2: Longitudinal momentum distribution of hadrons in 
a jet originating from a 20 GeV quark in vacuum and in 
medium. The medium results are given for two different val- 
ues of line-integrated virtuality AQ 2 for scenarios (A,B,C) 
describing three different characteristic paths through the 
medium (see text for details). 



energy from hard partons to soft radiation increases with 
the strength of the medium transport coefficient. 

One notes that for small AQj ot , there is good scaling be- 
tween all scenarios, and only for larger values of AQ\ ot 
the precise pattern at which time the virtuality is trans- 
ferred to the shower starts to be significant and longer 
paths lead to more suppression of the high qx tail. How- 
ever, since partons close to the edge of the medium prop- 
agating outward traverse only low density matter, they 
probe a region in which the scaling holds. Thus, scenario 
(B) can effectively be assumed to be a fair approximation 
for all paths if AQ\ ot is used to characterize the path. 
This observation simplifies the subsequent averaging of 
paths through the hydrodynamical medium considerably, 
as it reduces the problem to determining AQ\ ot on line 
integrals through the medium. 



E. The Hump-backed plateau 

In order to focus more on the hadron production at low 
p z , we introduce the variable £ = ln(l/.x) where x = 
p/Ej et is the fraction of the jet momentum carried by a 
particular hadron and Ej et is the total energy of the jet. 
The inclusive distribution dN/d^, the so-called Hump- 
backed plateau, is an important feature of QCD radiation 
[3r3 . [37} and is in vacuum dominated by color coherence 
physics. 

Fig. [3] shows the calculated distribution in vacuum and 
its distortion by the medium effect. Qualitatively, the 
results are very similar to those observed in (2(| and [2l[ 
for schematic radiative energy loss, i.e. a depletion of the 
distribution at low £ and a strong enhancement of the 
hump. 



I 4 



• vacuum 

• AQ 2 = 3GeV 2 

• AQ 2 = 10 GeV 2 



- 

* • 
* * 

* 1 i_ 



FIG. 3: The Hump-backed plateau distribution dN/d£ for a 
jet originating from a 20 GeV quark shown both in vacuum 
and medium modified at two different values of line-integrated 
virtuality AQ 2 . 



Transverse jet momentum spectra and angular 
distribution 



Since the parton shower picks up virtuality from the 
medium, we may expect some broadening of the trans- 
verse distribution of hadrons in the jet due to the medium 
effect. We show the transverse momentum qx distribu- 
tion (where transverse in this section denotes the direc- 
tion transverse to the hard parton initiating the shower, 
not as in a description of the whole p-p or Au-Au col- 
lision, the direction transverse to the beam axis) of 
hadrons in the shower in Fig. [¥] 
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FIG. 4: Momentum distribution of hadrons in the jet trans- 
verse to the direction of the shower initiating parton for a 20 
GeV quark, both in vacuum and in medium for two different 
values of the line- integrated virtualtiy AQ 2 . 

Although the effects on the transverse momentum spec- 
trum are not large, there is a clear trend to a rise of 
multiplicity at low qx and some depletion in the high qj- 
tail through the medium effect. One may also note that 
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the integral f dqxqrdN / dqx increases by ~ 30% from 
vacuum to the AQ| ot = 10 GeV 2 medium modified jet. 
The effect of the medium on the transverse structure of 
the shower is more clearly seen when looking at the an- 
gular distribution of hadrons in the jet. The distribution 
dN/da where a is the angle between hadron and jet axis 
is shown in Fig. [5] where a cut in momentum of 2 GeV 
has been applied to focus on hadrons which would appear 
above the soft background of a heavy-ion collision. 



mentation function in the hydrodynamically evolving 
medium. The nuclear suppression factor is defined as 



Raa(Pt,v) = 



dN h AA jdP T dy 
T AA {0)daPP/dP T dy 



(21) 



where TaaO 3 ) is the standard nuclear overlap function. 
We can compute it by forming the ratio 
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FIG. 5: Angular distribution of hadrons with respect to the 
jet axis for a shower initiated by a 20 GeV quark, both in 
vacuum and in medium for two different values of the line- 
integrated virtualtiy AQ 2 . 

The angular broadening of the distribution by the 
medium is clearly visible (note that the dip at a = 
is caused by the Jacobian in a polar coordinate system 
and is not an indication of a double-hump angular struc- 
ture). 



Raa(Pt,v) 



, AA-^h + X rpp- 
a0 medium l a " 



■h+X 



dP^dy 



dPldy 



(22) 



where da pp ^ h+x / dP^dy follows from Eq. Q when 
Df-th{z, (J if) is set to be the vacuum fragmentation func- 

/dP^dy is computed from the 



tion whereas da^+ x 



same equation with Df^h(z, [A) replace by the suitably 
averaged MMFF (D m m {z , fJ%)) t aa ■ This averaging has 
to be done over all possible paths of partons through the 
medium. We make the simplifying assumption that in- 
stead of computing the detailed paths of all partons in a 
shower through the medium, we only sample the medium 
along the eikonal path of the shower originator and con- 
sider this the medium seen by all partons in the shower. 
In essence, this neglects the transverse spread of paths in 
a shower. 

The probability density P(xo,yo) for finding a hard ver- 
tex at the transverse position ro = (£o>2/o) & n d impact 
parameter b is given by the product of the nuclear profile 
functions as 



r A (ro + b/2)r A (r -b/2) 
T A A(h) 



(23) 



IV. COMPUTATION OF Raa 

In this section, we show how to compute the nuclear 
suppression factor Raa using the medium-modified frag- 



where the thickness function is given in terms of 
Woods-Saxon the nuclear density p^(r, z) as Ta{v) = 
J dzp A {r,z). The MMFF must then be averaged over 
this quantity and all possible directions <f> partons could 
travel from a vertex as 



J 



(D M m{z^ 2 ))t a/ 



2tt jo 



dx / dy P(xo,yo)D M M(z,p 2 X)- 



(24) 



r 



Using the approximate scaling relation described in sec- 
tion [niDl the med ium modified fragmentation function 
Dmm(z,H 2 X) f° r a path £ can be found by computing 
the line intergal 



aqL = 



dCq(C) 



(25) 



through the hydrodynamical medium where <f(£) is given 
by Eq. (fl8|) . The MC shower code is then used to com- 
pute D MM (z, fi 2 , for each value of AQf ot . 

Note that this point that there is in principle a concep- 
tual difference between using a fragmentation function 
extracted from a shower simulation (be it in medium or 
vacuum) and a phenomenological fragmentation function 
such as the KKP |I| or AKK [H set in Eq. ©. This has 



8 



to do with the fact that in the the usual fragmentation 
functions \i t represents the relevant hadronic momentum 
scale. On the other hand, one can extract fragmentation 
functions from a shower evolution only for a given par- 
tonic momentum scale. Thus, the scale evolution of the 
fragmentation function must be accounted for in a differ- 
ent way. At present, we use a MMFF for a fixed partonic 
scale of 20 GeV, hence the description of high Pt hadron 
momentum spectra with the shower generator extracted 
fragmentation function is not as good as using the phe- 
nomenological functions. However, here we are chiefly 
interested in the ratio of modified over unmodified spec- 
tra, and in this ratio the uncertainties due to the scale 
evolution largely cancel out. 



V. COMPARISON WITH DATA 



MMFF model finds a drop. In particular, the region 
between 10 and 14 GeV is not well described. This is 
quite unlike models based on energy loss of the leading 
parton where a rise of Raa with pt is found, cf. e.g. 

[12]. On the other hand, a similar drop of Raa has 
been found in the schcmatical shower evolution model in 




It seems clear that there is a direct link between the 
characteristic distortion of the Hump-backed plateau at 
larger £ by medium effects and the rise of Raa towards 
lower pt — in both cases, the effects of additional parton 
production at lower momenta by medium induced radi- 
ation become visible. It seems at this point that this is 
a generic, unavoidable feature of the model formulation. 
Nevertheless, let us revisit the assumptions of the model 
and discuss their validity. 



At this point, our the model is completely determined 
except for the parameter K in Eq. (fT8|) which links the 
energy density in the medium with the transport coeffi- 
cient. We adjust K such that the best fit to the data is 
achieved. In Fig. ® we show Raa for different values of 
K. 



• PHENIX data 

— K= 1.0 

— K= 1.25 

— K= 1.5 
- K = 1.75 




) ' 1 ' 1 ' ' — 

5 10 15 

p T [GeV/c] 



FIG. 6: The nuclear suppression factor Raa for central AuAu 
collisions at 200 AGeV computed for different values of the pa- 
rameter K specifying the transport coefficient of the medium 
in comparison with the PHENIX data 0|. 

The data seem to favour K = 1.5. This implies that the 
medium exhibits a 50% stronger effect than expected in 
an ideal QGP [4fJ. qo, i.e. the highest transport coef- 
ficient reached in the evolution in the medium center at 
thermalization time is found to be 7.8 GeV 2 /fm. Note 
that this is well below the values of K (and q) extracted 
from computations which consider energy loss from the 
leading parton in the ASW formulation of energy 
loss, cf. where K values between 2 and 5 are found 

However, it is evident that the overall trend of the data is 
not well described by the model. While the data indicate 
a small rise of Raa from px = 5 GeV and above, the 



VI. DISCUSSION OF MODEL ASSUMPTIONS 

While a detailed discussion of all model uncertainties is 
beyond the scope of this paper and will be postponed 
to a future publication, we will at this point provide a 
list of possible uncertainties and estimate their effect if 
possible. 

First, there are uncertainties in the calculation of the 
hard process. We do this in LO pQCD, and one may 
wonder about NLO corrections. However, in (42j is has 
been shown that LO pQCD supplemented with a K fac- 
tor to account for higher order effects is in fact a good 
description of hard hadron production in p-p collisions. 
Another approximation we have made is to neglect the 
scale dependence of the fragmentation function and in- 
stead to compute the fragmentation function at a single 
partonic scale of 20 GeV. While there is no conceptual 
problem with a MC computation of the dependence of 
the vacuum or medium-modified fragmentation on [z, fi) 
and in the latter case also AQ^ ot , mapping out a depen- 
dence on three parameters is very demanding in terms of 
MC simulation time and will be done in a future investi- 
gation. As shown in Fig.[7]the effects are not substantial. 
A second class of uncertainties has to do with the way 
the shower algorithm is implemented and affects both 
vacuum and medium results. In essence, these are en- 
coded in the control parameters of the PYTHIA shower 
algorithm. For example, the scale Q^in down to which 
partons are evolved is a parameter, and the results show 
some dependence to it. Other uncertainties have to do 
with what scale is used to determine the running of a s 
in the shower. In general, these effects are known to be 
small. 

A third class of uncertainties has to do with the spe- 
cific implementation of medium effects. Chiefly, one may 
wonder about the effect of replacing a probabilistic for- 
mation time with its average value in Eq. (| L6[) or about 
replacing a probabilistic transfer of virtuality associated 
with scattering processes with an average growth of vir- 
tuality as realized in Eq. (|17p . The latter assumption 
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FIG. 7: Dependence of the fragmentation function, both in 
vacuum and in medium for a line-integrated virtuality AQ 2 = 
10 GeV 2 for two different values of the partonic momentum 
scale p,. 



has some impact on the treatment of QCD coherence ef- 
fects, as the angular ordering Eq. (fTo]) would be destroyed 
by a scattering process. Thus, a detailed formulation in 
terms of individual scatterings with medium constituents 
would include the information when angular ordering is 
effective and when not. This information is lost when 
computing based on an average virtuality transfer. On 
the other hand, Landau-Pomcranchuk-Migdal (LPM) in- 
terference [43[ would lead to a suppression of radiation as 
compared to the vacuum situation due to frequent inter- 
actions with the medium. While answering the question 
about the effect of replacing a probabilistic formulation 
with averages requires a restructured algorithm which is 
beyond the scope of this paper, the effect of a detailed 
treatment of the angular ordering can be estimated com- 
paratively easy by computing in the two limits of re- 
taining or dropping angular ordering completely in the 
medium. As evident from Fig. [8l the effects are in prac- 
tice relatively small and moreover do not alter the shape 
of the distribution significantly, hence they can largely 
be absorbed into a rescaling of AQ^ ot or alternatively K. 
One may also keep in mind that using the Lund string 
fragmentation mechanism to hadronize in a heavy-ion en- 
vironment is not straightforward, especially at low z and 
for heavy hadrons. If one estimates the hadronization 
time for a hadron with mass rrih and energy Eh as 



Eh /m 2 h 



(26) 



then a 10 GeV pion will be formed several tens of fm away 
from the medium, but a 5 GeV proton has a formation 
time of about 1.2 fm, clearly not enough to assume the 
process happens out of the medium or even out of the 
partonic region. Thus, there is some reason to suspect 
that the treatment of low z hadrons will miss out part of 
the relevant physics. 

Finally, one may question the validity of the scaling as- 
sumption made in section IIII Dl which allows in essence 
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FIG. 8: Longitudinal momentum distribution of hadrons in a 
jet originating from a 20 GeV quark in vacuum and for two 
different values of the line-integrated virtuality AQ 2 . Shown 
is the effect of retaining angular ordering in the branching (as 
in vacuum) or dropping it under the assumption that frequent 
medium interactions destroy coherence. 



to reduce the information of q(£) specified on the whole 
path to a simple line integral. While the scaling is not 
exact, a detailed comparison of Eq. (fT7| for some selected 
explicit forms of paths through the hydro medium with 
the scaled version shows that differences in the resulting 
MMFF are at most on the level of 20% and moreover can 
again largely be absorbed into a rescaling of AQ 2 ot as the 
shape of the MMFF is not changed. 



VII. CONCLUSIONS 

There is evidence that model uncertainties which can be 
estimated maninly affect the precise determination of K, 
but not the more striking fact that model does not cap- 
ture the rising trend of Raa with px apparent in the data. 
While there is still the possibility that one of the effects 
which have been neglected in modelling so far may cor- 
rect this disagreement, one may note that this does not 
seem likely as there is within the model a rather generic 
physics mechanism which produces this particular shape, 
i.e. the energy transport from high z partons into radi- 
ated low z multiplicity, and it is this additional multi- 
plicity which leads to the rise of Raa for low momenta. 
Supporting evidence for this is found in the fact that the 
calculation in [2(| which is based on a completely differ- 
ent algorithm but models the same physics mechanism 
results nevertheless in a very similar shape of Raa when 
folded with the pQCD parton spectrum. 
At this point, we note that there's a second failure of 
the model, also connected with subleading hadrons in 
the shower. The angular distribution of hadrons shown 
in Fig. [5] shows, even for a jet emerging from the center 
of the medium, only moderate broadening, but no sign 
of an energy momentum flow to large angles or the for- 
mation of a dip at zero angles as would be required to 
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account for experimental data on two and three particle 
correlations |4{| H|| [13] • In contrast, models which 
assume that the energy radiated from the leading parton 
is not developing as a shower but rather excites a hydro- 
dynamical Shockwave in the medium [48], [4t| [5(J HH can 
account for the observed structures in the data. 
On the other hand, radiative energy loss considered for 
the leading parton only is able to account well for the 
data even in the nontrivial case of back-to-back correla- 
tions [12) HH, E3- Furthermore, conceptually branchings 
at large Q 2 should take place regardless of the presence or 
absence of a medium, as they on average develop before 
a thermalized medium can be assumed to exist. 
These observations lead to the idea that while especially 
the initial development of a high pt process is described 
by the shower evolution, at some lower scale a pQCD 
shower ceases to be an adequate description for soft ra- 
diated partons, but instead some form of a coupling to 
hydrodynamical modes of the medium (cf. e.g.[52f) takes 



place. The interesting question (especially in the light of 
the kinematic range of hard probes accessible at LHC) is 
what part of the shower (starting from the leading par- 
ton) is still correctly described by the in medium shower 
evolution equations and under what conditions this evo- 
lutions should be replaced by different physics. The an- 
swer to these questions at RHIC energies may be found in 
the application of the formalism presented here to hard 
multiparticle correlations in which a systematic access 
to the high Pt part of the shower is accessible. This, 
however, will be the subject of a future investigation. 
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